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Abstract. Recent observations of the Galactic center revealed a nuclear disk of young OB 
, ^ . stars near the massive black hole (MBH), in addition to many similar outlying stars with 

higher eccentricities and/or high inclinations relative to the disk (some of them possibly 
belonging to a second disk). In addition, observations show the existence of young B stars 
(the 'S-cluster') in an isotropic distribution in the close vicinity of the MBH (< 0.04 pc). We 
use extended N-body simulations to probe the dynamical evolution of these two populations. 
We show that the stellar disk could have evolved to its currently observed state from a thin 
disk of stars formed in a gaseous disk, and that the dominant component in its evolution 
is the interaction with stars in the cusp around the MBH. We also show that the currently 
observed distribution of the S-stars could be consistent with a capture origin through 3- 
body binary-MBH interactions. In this scenario the stars are captured at highly eccentric 
orbits, but scattering by stellar black holes could change their eccentricity distribution to be 
C^ . consistent with current observations. 
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1. Introduction contain many O-stars residing in a stellar 

C*~) ■ disk moving clockwis e in the gravitational 

High resolution observations have revealed the poten ti al of the M BH ( Levin & Beloborodov, 

existence of many young OB stars in the galac- 2003t Genzeletal.. 2001 iLuetal.. 2006; 

tic center (GC). Accurate measurements of the iPaumard et all |2006|; |Tanneret alj,|2006|). The 

orbital paramters of these stars give strong evi- orbits of the stars in this disk have average ec- 

dence for the existence of a massive black hole centricity of ~ 0.35 and the opening of the 

> ! (MBH) which govern the dynam i cs in t he GC disk is h l R ~ 0.1, where h is the disk height 

dSchodel et all 120021: iGhez et aUl2003l) . Most and R is its r adius [Paumard et al. (2006) and 

of the young stars are observed in the cen- |Tanner et al.| (|2006|) suggested the existence of 

tral 0.5 pc around the MBH. The young stars another stellar disk rotating counter clockwise 

population in the inner 0.04 pc (the 'S-stars' and almost perpendicular to the CW disk. This 

or the ' S-cluster' ) contain only young B-stars, disk is currently debated as many of the young 

in apparently isotropic distribution around the stars are have intermediate inclinations, and are 

MBH, wit h relatively high eccentricities (0.3 ;£ possibly just outli ers that do not f orm a coher- 

e < .95) dGhez et alll2003llEisenhauer et all ent disk structure (|Luetal.|,[200g). 



2005). The young stars outside this region 
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Here we briefly report on our study of the 
dynamical evolution of the young stars in the 
GC, both in the stellar disk and in the S -cluster. 
We use extensive N-body simulations with re- 
alistic number of stars (10 3 - 10 5 ) using the 
gravitySimulator, a 32-node cluster at the 
Rochester Institute of Technology that incor- 
porates GRAPE accelerator boards in each of 
the nodes dHarfst et all 120071) . Thus we are 
able to probe the dynamics of the stars near the 
MBH and their stellar environment. We study 
two basic issues: (1) the long term evolution 
of the S-stars up to their lifetime of a few 10 7 
yrs, including their dynamical interaction with 
stars in the vicinity of the MBH; (2) The evo- 
lution of a realistic stellar disk, taking into ac- 
count both the effects of non-equal mass stars, 
as studied earlier, and more importantly the ef- 
fect of the interactions of disk stars with the 
stellar cusp around the MBH. As we show, the 
latter component proves to be more important 
than the other components discussed in previ- 
ous studies. A detailed report of our complete 
set of out simulations (not shown here), in ad- 
dition to analytic calculations will be presented 
in upcoming papers (Perets et al., in prepara- 
tion) 

2. Formation and/or migration origin 

Analytic calculations and simulations have 
shown that a young stars could have formed 
and grown over short times of thousands to 
millions of year s in a gaseous disk around 
the MBH (e .g. dNayakshin & CuadraL 120051: 
iLevinl 120071) ). Such stars could then form the 
stellar disk currently observed in the GC. It 
was suggested that the 'S-stars' with their very 
different properties migrated from the stellar 
disk t hrough a pla netary migration like pro- 
cess dLevinl 120071) . This interesting possibil- 
ity has not yet been studied quantitatively, but 
would suggest that the migrating stars should 
have relatively low eccentricities. Another pos- 
sibility is that these stars have a different ori- 
gin, possibly from the disruption of young bi- 
naries and the fo llowing capture of one of 
their components dGould & Ouillenl 120031) . It 
was recently shown that such a scenario could 
be consistent with the current knowledge re- 



garding the number of the observed 'S-stars' 
dPerets et all 120071) . The initial eccentricity of 
the captured stars should then be very high 
(> 0.96) in this scenario. We note that other 
scenarios were suggested for the origin of 
the young stars in the GC, but seem to be 
disfavored by cu rrent o bservations (see e .g. 
lAlexan der (2005) and Pa umardet all (12006b '). 



3. The S-stars 

Only one study has been published of the 
dynamical evolution of the S-stars since 
their capture / formation, which explored the 
possi ble role on I MBH on their evolution 
dMikkola & MerrittL 120081) . Here we present 
the results of N-body simulations dedicated to 
study such evolution. In order to do so we mod- 
eled a small isotropic cusp of 1200 stars, with 
masses of 3M (200 stars) and 10M o (1000) 
around a MBH of 3.6 x 10 6 M Q . We used a 
power law radial distribution of r~ a extend- 
ing from 0.001 - 0.05 pc near the MBH, with 
a — 2 for the more massive stars and a - 
1.5 for the lower mass stars. The more mas- 
sive stars correspond to the many stellar black 
holes (SBHs) thought to exist in this region, 
whereas the lower mass stars correspond to the 
S-stars in the same region. Since some of the S- 
stars may have higher masses of ~ 10 M , the 
higher mass stars in the simulation could also 
be treated as S-stars. We did not see any major 
differences in the evolution of the more mas- 
sive and the less massive stars, and we discuss 
the evolution of both together. 

We studied two evolutionary scenarios for 
the S-stars. In the first we assumed that the S- 
stars were captured thr ough a binary disruption 
scenario by the MB H dGould & Ouillenl 120031: 
iPerets et all 120071) and therefore have initially 
highly eccentric orbit (> 0.96) and they evolve 
for few 10 7 yrs. In the second scenario we as- 
sumed the S-stars formed in a gaseous disk and 
migrated to their current position, and there- 
fore have low eccentricities (< 0.3) and they 
evolved for 5 Myrs (the lifetime of the ob- 
served stellar disk. In order to check both sce- 
narios we followed the evolution of those stars 
in our simulation with highly eccentric initial 
orbits (the first scenario) and those with low 
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eccentricities (the second scenario) for the ap- 
propriate time scales. In Fig. (1) we show the 
final eccentricity distribution of the S-stars in 
both scenarios, as compared to the the orbits of 
the observed S-stars (taken from Gillessen et 
al.). These results suggest that, given the small 
statistics (16 S-stars with known orbits), the 
first scenario is much favored since it could 
explain the currently observed orbits of the S- 
stars, i.e. stars on highly eccentric orbits could 
be scattered by other stars or SBHs to smaller, 
and even much smaller eccentricities during 
their lifetimes. The second scenario, however, 
seems to be excluded (for the given assump- 
tions), since it has major difficulties in explain- 
ing the large number of eccentric orbits in the 
S-stars observations vs. the bias towards low 
eccentricity orbits seen in the N-body simula- 
tions simulations. This is clearly seen both af- 
ter 5 Myrs of evolution, and even after longer 
evolution, if these stars formed in an earlier 
epoch of star formation in a disk (not currently 
observed) 20 Myrs ago. 

4. The disk stars 



| Simulations: 20 Myrs 
I Observations 



lAlexander eTaTl d2007l) and ICuadra et all 
(2008) explored the dynamical evolution 
of a single stellar disk using small N-body 
simulations (~ 100 stars), where they studied 
the effects of massive stars in the disk (the 
mass function), and the structure of t h e disk 
(eccentric vs. circular). ICuadra et al.l J2008I) 
also studied t he role of wide binaries following 
iPerets et all d2008l) who suggested binaries 
could have an important role in the evolution 
of the disk, some what similar t o their role 
in stellar clusters (Heggig, 1 19751) and in the 
ejection of OB runaway stars. These studies 
showed that although the different components 
contribute to the disk evolution, it is difficult 
to explain the current eccentricities of the 
observed stars and the thickness of the disk, 
with only these components. 

We studied a single disk of ~ 5000 M , 
composed of either 5000 equal mass stars, or 
~ 2500 stars with a Salpeter mass function 
between 0.6 - 50 M . The initial conditions 
are of a thin disk (H/R ~ 0.01) with a surface 
density of r~ 2 , with all stars on circular 




^^Simulations: 20 Myrs 
| | Simulations: 5 My: 
^H Observations 




Fig. 1. The eccentricities of observed and simulated 
S-stars from the binary disruption scenario (after 5 
Myrs; upper figure) and from the disk migration sce- 
nario (after 5 and 20 Myrs; lower figure). 



orbits (e < 0.01). In addition we studied the 
evolution of such disks both with and without 
an isotropic stellar cusp component around the 
MBH in which the stellar disk is embedded 
in. The region of the GC disk is thought to 
contain a few 10 5 up to 10 6 stars; simulating 
such a larger number of stars in orbits close 
to a MBH, even with a GRAPE cluster is 
currently difficult. However, the dynamics 
and relaxation processes close to the MBH 
are dominated mostly by the much smaller 
number of SBHs in this region (a few 10 3 up 
to 10 4 SBH s are though to exist in the GC 
disk region :Miralda-Escude & Gould, 2000; 
Hopman & Alexander! 120061 : iFreitag et all 
2006). Simulating only this SBHs component 
could therefore be more efficient in running 
time and at the same time capture most of the 
important relaxation processes effecting the 
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the dynamics of the disk. In our simulations 
we put 1.6 x 10 4 SBHs (of 10 M each) with 
an isotropic power law distribution (n oc r~ 2 ), 
between 0.01 - 0.8 pc. 

The evolution of the mean eccentricity of 
the disk stars is shown in Fig. (2), both for 
low mass stars (< 15 M ) and high mass stars 
(M > 15 M Q , i.e. such as the observed disk 
stars in the GC). The evolution of a stellar disk 
embedded in a cusp of SBHs, is compared with 
that of an isolated stellar disk with a Salpeter 
mass function between 0.6 - 80 M (i.e. higher 
mass cutoff than used in the disk+cusp simu- 
lation, to allow for the disk heating by more 
massi ve stars, as discussed in lAlexander et all 
120071) . The results show that the SBHs cusp 
component has an important contribution to 
the disk evolution. Although one disk has a 
lower mass cutoff than the other, it is heated 
much more rapidly, due to the contribution 
of the cusp stars. More importantly, the high 
mass stars corresponding to the stars that are 
actually observed in the GC disk, have rela- 
tively low eccentricities in the isolated disk, 
and would present difficulty to our under- 
sta nding of the di s k evolution, as discussed 
by lAlexander et all l2007t ICuadra et all l2008l 
Adding the cusp component solves the prob- 
lem, as even the eccentricities of the higher 
mass stars are high in this case and compara- 
ble with the observed eccentricities of the disk 
stars in the GC. We note that these simulations 
did not take into account the contribution of 
low mass disk stars (i.e. not the SBHs in the 
cusp) which will further accelerate the heating 
of the stellar disk. 



5. Summary 

The dynamical evolution of the young stars in 
the GC both in the stellar disk(s) and in the 
inner S-cluster is not yet understood. We used 
N-body simulations to study the dynamics and 
origin of these stars. We found that the S-stars 
close to the MBH in the GC could be stars 
that were captured following a binary disrup- 
tion by the MBH, and later on dynamically 
evolved due to scattering by other stars, or stel- 
lar black holes, to obtain their currently ob- 
served orbits. We also show the the young stel- 
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Fig. 2. The evolution of the mean eccentricity of the 
disk stars (both low mass, M < 15 M Q and high 
mass stars, M > 15 M ) with and without interac- 
tions with the stellar cusp around the MBH. 



lar disk could have formed as a cold (thin) cir- 
cular disk and evolve to its currently observed 
thick (hot) disk, mostly due to scattering by 
cusp stars, whereas self relaxation of the disk 
plays a more minor role, especially in regard 
to the more massive stars seen in observations. 
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